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ABSTRACT 


The validity O£ Conventional skin and stringer type 
analyses is verified experimentally for a thick milled skin 
atecrart Structure. The accuracy of solution is evaluated 
by comparing the strain energies in a mathematical model to 
that obtained experimentally from a wiffig Paneer eamnigh 
performance aircraft. 

Results indicate that th®-conwerttional methods of a 
skin-stringer/matrix force approach are valid for a thick 
milled skin aircraft structure of this type provided the 
panel loading is neither severe nor concentrated. Addi- 
tionally the grid elements chosen must be small enough to 


preclude any gross averaging errors. 
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CHAPTER TI 


INTRODUCTION 


Aircraft structures in general, and wings in parti- 
cular, have been constructed using the conventional stringer 
and sheet methods. This form of construction, by nature not 
a simple structure, introduced the problem of how to accu- 
rately predict the load distribution and especially the 
stress concentrations that could be developed. The oldest 
representation scheme for the discrete element analysis of 
the problem evolved as the axial force member-shear panel 
idealization that is described in Reference 1. This method 
1s particularly aceestable to the aircraft industry as it 
readily lends itself to various matrix techniques that are 
easily adaptable to a digital computer for solution. Gen- 
erally speaking, this method of solution makes no provision 
for accounting for Poisson's ratio and sweep effects. 

With the advent of the very heavy milled skins of 
modern high performance aircraft, one immediately questions 
the validity of the conventional form of solution. In order 
to properly account for bending-torsion interaction, chord- 
wise curvatures, sweep coupling, shear deflections, and the 
increasing use of large cutouts, several alternative methods 
of analysis have been developed. 

The most widely used solution technique is outlined by 
J. H. Argyris and S. Kelsey in Reference 2. The aircraft 


structure is molded into a large assemblage of elastic ele- 


ments, which permits matrix formulation of the solution in 
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terms of the different energy theorems. The fundamental 
principles originated by Maxwell and Mohr of virtual forces 
and virtual displacements are used in the energy method of 
Seelucion. A practical application of this technigue is 
included in References 3 and 4. 

The objective of this study is to determine the valid- 
ity of this lumped parameter form of a theoretical analysis 
for thick skinned structures. Because of the general energy 
approach to the solution, the method chosen to determine the 
accuracy Of the results was to compare the total internal 
strain energy of an idealized panel with 8 OL dh “actual 
thick milled skin aircraft structure of the same dimensions. 
mnie Struceurc wOUld wave veo De fhmick Gnough tombe aaleeco 
SUpPpOrt” ay Dehd Ing erorsiGnm Interaction that may Gxist im 
addition to the inplane loadings. 

The objective of this report was to be achieved by 
comparing the theoretical strain energy of a structure with 
that determined empirically from strain gage readings. The 
structural specimen to be analyzed was a single wing panel 
taken from the port wing upper skin of a Mach 2.5 all- 
weather fighter project which was canceled by the U. S. 
Navy before the structural test program was completed. 

The necessary strain gage rosettes and accompanying 
instrumentation were installed, and a comparison of the 
theoretical and empirical analyses was made. 

A literature search did not reveal any reference that 


made a comparison of the energy level in an idealized 
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structure to that of an actual wing panel element. The 
importance of a comparison of this type is apparent when one 
realizes that the integrated behavior of a structure is 
dependent upon the energy accounting in the system. 

This study was completed during the 1967-1968 academic 
year at the Naval Postgraduate School, Monterey, California. 
Acknowledgement is gratefully made to Professor C. H. Kahr 
of the Aeronautics Department for his guidance and consulta- 
tion as thesis advisor; R. A. Besel and T. B. Dunton, 
laboratory technicians; and CAPT D. B. Messerschmidt, USMC 
who, while working on a similar thesis project, aided in 


Enis study. 


1. 


CHAP TE Raa 


THEORETICAL ANALYSIS 


The structure chosen for study was the wing center 
Secttonesot an FSU-3 type aircraft. The particular panel 
analyzed was the port section bounded longitudinally by the 
pivot and intermediate ribs and laterally by the aft- 
intermediate and rear-intermediate beams. The wing center 
section and specific panel are graphically depicted in 
Figure l. 

The typical industrial approach to poeey 21d SMe iG= 
ture of this nature would be to determine what valid sim- 
plifying assumptions could be made in order to create a 
mathematically simple idealized structure. The resulting 
model would consist of axial load carrying bara and plane 
shear carrying constant thickness panels. The idealized 
model was developed and is thoroughly discussed in 
References 3 and 4. 

The load application was accomplished by applying a 
torsional moment at the wing fold rib, about the load 
reference axis (CIB). Considering the port wing inboard of 
the intermediate rib, this couple must be distributed, 
using elementary torque theory, to the equivalent loads on 
a streamwise cross-section (intermediate rib) at station 
y,, = 81.98 (Figure 2). The loading was then divided among 
the six cells of the wing. The shear flow for each cell 


was determined in accordance with Reference 1 by equating 
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FIG. | 
SKETCH OF F8U-3 WING CENTER SECTION 


FIG. 2 
CENTER SECTION — F8U-3, PORT WING 
(Main structural elements) 
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Wing station, X,,, in. 


the angle of twist of one cell with the remaining five and 


equating torsional moments as follows: 
Tihcanaeuh 
= (1) 
, 4 


The results and calculations of the shear flows were 
Calculated in Reference 3 and are included in this report 
as Appendix A with inappropriate sections deleted. Using 
the planform dimensions of Figure 2, the shear flow for each 
cell was determined as oriented perpendicular to the CIB. 

Figure 3 shows the idealized panel section of the wing 
with pertinent dimensions obtained from Reference 5. The 
internal strain energy of this skewed shear panel was 
calculated in accordance with S. J. Garvey (Reference 6) as 
follows: 

"Since by similar triangles, the perpendicular distance 
of any point in a trapezoidal panel to a parallel edge is 
proportional to the depth of the panel at that point 
measured parallel with a side, the well-known result follows 
that the "Shear flow” in a trapezoidal panel is inversely 
proportional to the square of the depth. "© 

To find the strain energy, using these geometric sim- 
plifications of a trapezoid, one must solve the following 


equation: 


U= ste f 7 sm y « Volo — 


i 
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where: 


v= cot @ + cotg@ cot & + cot SC (3) 


Table I lists the appropriate dimensions and constants 
needed for the theoretical analysis. The volume of the 


plate was determined for the idealized panel from: 
V = as “Ste ) b (tr 7 t. | (4) 
= 


This approach of utiliZing the average thickness will ob- 
viously limit the accuracy of the resulting theoretical 
strain energy; however, the objective of this part of the 
analysis is to obtain the strain energy by tne same conven- 
tional methods currently being used by industry. 

As previously discussed, the plate edge loadings were 
determined from a linear interpolation using the planform 
dimensions. These shear panel edge loadings were distri- 
buted around the perimeter of the panel, in accordance with 
Reference 4, by first assigning a unit shear force to the 
Outer edge of each panel. Considering a typical panel as 
shown in Figure 4 and taking moments about appropriate 
points to ensure the maintenance of equilibrium, the follow- 


ing ratios of edge forces were derived: 


. 5 
ap = a. tM (5) 
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TABLE I 


PANEL PARAMETERS 


O52 7079—=76 Al 


3.9 x 106 7079-T6 Al 


10.3 x 10° i. 7079-T6 Al 


L232 Boi. 


Lome Rebs 


365 605 l ikeces 


Ref. 
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4 Aa (7) 
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The Matrix Force Method is described in Reference 11 
and was used to determine the end edge loadings on the 
panels. These loads were obtained from previous calcula- 
tions in Reference 3. The end edge load for the particular 
panel studied in this report (Figure 5) was then scaled to 
the edge utilized in the strain energy calculation of 
Equation (2). This scaled load was then divided by the 
length of the same edge. The result was the desired edge 
shear flow (q, = 325.976 pounds per inch). 

Substitution of the appropriate values from Table I and 
the computed shear flow into Equations (2), (3), and (4) 
will yield the total shear panel internal strain energy. 
The resulting value was 33.280 inch pounds. 

The idealized model, being a pure shear panel, was in- 
capable of supporting any axial loadings. For this reason 
the effective cross sectional area of the plate was divided 
evenly and distributed to the adjacent flanges (Figures 3 
and 5). The axial loads were then assumed to be applied 
directly to the flanges. 

In accordance with Reference 9, sets of additional 
axial forces (dashed forces in Figures 3 and 5) were applied 


to the ends of the flanges. The forces are necessary to 
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account for the additional strain energy stored by the axial 
components of shear flows caused by the sweep effect of the 
non-rectangular panels. It is important to realize that 
these sweep coupling loads are internal forces and would not 
Demimeiuded in the equilibrium equations) of |\the structure. 

From an energy standpoint, these dashed forces account 
for the additional strain energy stored in the panel. The 
energy is accounted for in this manner because the longi- 
tudinal flanges contain the effective area of the panel. 

The contribution of a uniform bar under linearly vary- 
ig axial *orees to ‘thes total iIntemnal straingeenersy is 


given in Reference 9 as, 


£ A 
U = Be |(P + oe a poy x) d x (8) 


Where @Eee oid P; are the axial forces applied to the ends 


of the bar. Integration of Equation (8) yields: 


U= che (Po + P+ PB) 9 
The flange loadings were determined as a percentage of the 
total load based on the ratio of the effective flange area 
of the panel to the total flange area. 

Substitution of appropriate panel parameters and pre- 
viously determined flange loadings into Equation (9) yields 
the additional strain energies required to complete the 
theoretical energy accounting. > tie distribution @fene 


energies is listed in Table If. 
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TABLE IL 


THEORETICAL STRAIN ENERGY DISTRIBUTION 


OF THE IDEALIZED MODEL 


Shear Panel 33.280 inch pounds 
Forward Flange 0.313 inch pounds 


Rear Flange 4.067 inch pounds 


Rib Flange 0.046 inch pounds 


TOTAL ENERGY 34.706 inch pounds 
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CHAPTER. 1 1. 


DISCUSSION QF THEORETICAL RESULTS 


The concept of the idealization, as discussed in 
Retemence 1, is simply to assign all the direct stress 
Carrying capacity of the structural material to the axial 
force members and the shear carrying capacity to the 
shear panel. 

The idealization yields a theoretical solution to the 
problem that utilizes several different averaging processes; 
however, the objective of this report is to analyze the 
accepted conventional methods currently in use in the 
aircratt industry. 

The distribution of the panel area to the adjacent 
flanges is complicated by the apparent double use of the 
area near the tip of the panel. Several wing analyses have 
described alternative methods of including the rib flange in 
ielgver choles leave. peractire. It is of interest to note that 
the theoretical result of the strain energy determined for 
the rib flange is insignificant when compared to the other 
components of the strain energy in the structure. Because 
of the small energy of the rib flange portion of the plate 
contribution, one could conclude that this flange and simi- 
lar pieces do not contribute significantly to the over-all 
wing behavior. Therefore, as discussed in Reference 12, 
analyses of wings with applied normal loadings could he 
simplified considerably by assuming zero flexibility for 


this flange. It should be noted that the application of 
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Cchordwise loads would require that these elements be in- 
cluded. This approach would appreciably decrease the number 
of redundants in a wing analysis. 

An appropriate study might be the determination of the 
contribution of these rib flanges to the integrated analysis 


of the wing. 


Zo 


CHAPTER IV 


EXPERIMENTAL ANALYSTS 


iieomies tO determinesthewactiualestr ain enezoweand 
compare it with the value predicted from theory, an experi- 
mental investigation was completed on a panel of an F8U-3 
wing acquired from the Ling-Temco-Vought Corporation. The 
bemel Section was a Eypical thick milled skin portion ct 
the wing and was tapered along the span as well as in the 
chordwise direction. 

The detailed explanation of the wing installation and 
loading is included in Reference 3; therefore, only the 
highlights of installation and any modifications incorpora- 
ted by the author will be included in this report. 

The wing was mounted inverted on a rigid support jig. 
It was fastened to the jig at the two fuselage pivot lugs 
and the two points on the main box front beam directly 
opposite the bumper points (Plate 1). Aluminum pads were 
installed under the jig to distribute the loading to the 
avo rcacOn zee lOOr. 

The entire wing was subjected to two pure torque loads 
applied at each wing fold rib about the CIB (Plate 2). The 
position of the elastic axis was obtained from manufac- 
turer's data and is shown in Figure 2. The CIB is approxi- 
mately parallel to the elastic axis in the outboard section 
Of the wing...It reflects the characteristic breakdovwuein 


the concept of a load axis toward the root of a swept wing. 


30 





oul 





ye 


The predetermined loading (8,000 pounds or 336,000 inch 
pounds of torqu¢) was measured with 10,000 pound dynamo- 
meters connected in series with eachehydmaulie loading 
cylinder. They were graduated in 10OO pound increments but 
could be read accurately to ten pounds. ‘The final load 
magnitudes were determined with hydraulic manifold pressure 
readings related to pressure-load calibration Serves fOr 
the cylinders found in Appendix B. The load was well 
within the elastic limit of the wing structure but still 
of sufficient magnitude to yield adequate strain levels 
throughout the structure. 

All outer surface strain gage installations on the test 
panel were installed by the author. The internal gages 
were installed by laboratory technicians. All strains were 
measured with SR-4 strain gage rosettes of the type AR-7-2 
manufactured by the Baldwin-Lima-Hamilton Corporation. 
Appendix C is a list of all 102 gages, their gage factors, 
resistances and co-ordinate locations on the panel. The 
rosettes on the outer surface were installed back to back 
with the inner gages in order to give an indication of any 
differential bending effects through the thickness of the 
panel that may have existed. Figures 6 and 7 graphically 
depict the strain gage locations on the panel, and Plate 3 
shows the interior gage installation. A twenty channel 
Budd Company switching and balancing unit was connected to 
an externally powered Wheatstone bridge circuit. The bridge 


and all associated circuitry were an integral part of a 
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Budd/Datran Digital Strain Indicator (Model TC22) that read 


out strain directly in micro-inches per inch. Calibration 


of all equipment preceded each run. The electronic set-up 


is described in Reference 3 and is reproduced in this report 


as Appendix D. 


The procedures used in conducting all test runs were 


as follows: 


i. 


Because 


sun striking 


The digital indicator was calibrated before 
each run ensuring that sufficient time was 
taken between runs to preclude any errors 
being induced from hysteresis effects. 

The switching and balancing unit leads were 
plugged into twenty gage terminals at the 
junction panel (Plate 4). 

All strain gages and dynamometers were zeroed. 
The hydraulic torque loading was applied and 
adjusted to agree with all previous runs. 
Equal couple loadings applied to each wing 
tip were ensured by comparing the dynamometer 
readings. 

The strain readings were taken and recorded. 
The loads were removed. Hysteresis effects 
were removed by allowing sufficient time (8 to 


10 minutes) before commencing the next run. 


of ambient laboratory temperature changes, the 


the wing, and outside electrical interference, 


all test runs were conducted at night. To ensure steady 
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state conditions an adequate (approximately one half hour) 
SWESCPIE time was used for all associated equipment. Repeat- 
ability of strain gage readings was ensured by check runs 
being made periodically and compared with those previously 
obtained. The percentage difference in readings was minute 
and was considered to ave an insignificant effect on the 
analysis. 

As stated earlier, the objective of the experimental 
phase of this analysis was to correlate the theoretical 
results that were previously obtained; therefore, once the 
recorded strain data were obtained, the experimental total 
strain energy in the plate was to be determined. Refere..ou 
ence 7 discusses the calculation of strain energy ina 
structure in a general form that is applicable to the panel 
considered in this report. The general form of the total 


strain energy is given by, 


U= ff [se(@ + G+ )- G(aG+ Ga 
V 


(10) 
| 2 2 é 
+1 G)+ 26 (Tey + Tye + Tax) dy du d2 
In the wing panel of this report, where U3, Tuy2 
and (Pa are all nearly zero for plates/panels of this 


type, the stress condition indicated by Equation (10) 


reduces to the following: 


U=tf[[ee" ow, — ZU GGy) +o Tes) didy de = (11) 


6 


This strain energy formula of the stress state in 
Cartesian co-ordinates is the same as that related by 
A. G&G. M. Grzedzielski in Reference 8. Here the terms have 
the following meanings: Tk term represents the shear 
energy of the panel. The oq and <o reams reflect 
the inplane and bending/torsion energies. The term contain- 
ing Gx Gy introduces the Poisson's ratio coupling. 

Strain gage data were programmed into an IBM System/360 
Model 67 digital computer to obtain the principal stresses 
and directions at all rosette locations. A FORTRAN IV 
program "RECROSE"”" was written to convert the strain readings 
to the desired stresses. The discussion of the computer 
program and the results are included in this report as 
Appendix E. 

The test panel was then partitioned into fifteen sub- 
panel elements. The partitioning of the panel (Figure 8) 
would permit a summation of strain energies over the panel 
that would yield the total strain energy. This process is 
a close approximation to the integration required in 
Equation (11). The strain readings were averaged where more 
than one gage existed. Both the inside and outside rosettes 
on each sub-panel were included. The results of the sub- 
panel averaging are listed in Table III. The sub-panel 
thicknesses used to compute the volumes were the average 
sub-panel thicknesses obtained from Reference 10. A graph 
of the panel thickness versus y,, is included in this 


report as Figure 9. 
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TABLE It 


AVERAGE SUB-PANEL STRESSES 


1 = 905 25 
Z é -338.28 
3 -502.31 
= =O a Oe 
S =669.U6 
6 -691.19 
7 -898. 64 
8 -910.20 
9 =3/70. 08 
10 =—1l39.25 
lak : -1127.20 
12 =10642 52 
IL -1225.65 
14 =1 291296 


di -1285.93 
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The differences between the inner and outer surface 
stresses (Appendix E) indicate that a bending/torsion effect 
definitely does exist; however, assuming a linear variation 
of the stress through the thickness of the plate produces a 
strain energy that is less than one per cent different from 
a direct stress averaging technique. Therefore, the mid- 
plane stresses obtained by averaging the inner and outer 
surface stresses were used in the calculation of the experi- 
mental strain energies. 

The individual sub-panel strain energies were computed 
from Equation (11) using the required stresses listed in 
Appendix E. The sub-panel energies were then summed to 
obtain the total strain energy of the entire wing panel. 
The result of this summation was 33.714 inch pounds. 

The internal strain energy thus obtained was nearly 
that determined theoretically in Chapter II and was con- 
sidered a reasonably accurate evaluation of the actual 
strain energy that existed in the swept panel element. 

Figure 10 is a plot of the principal stresses obtained 
from "RECROSE" versus location on the panel. The purpose 
of this graph is to give the reader a feeling for the over- 
all panel loading under the applied loading of the whole 


wing. 
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CHARTER V 


DISCUSSION OF EXPERIMENTAL RESULTS 


The strain gages on the outer surface were mounted back 
to back with the interior rosettes. The purpose of this 
installation was to detect any bending present that would 
be reflected in a difference in inner and outer strain 
readings. Examination of Figure 10 and Table E4 in 
Appendix E indicates that bending definitely does exist. 
However, when a linear variation of stress through the 
thickness of the plate is assumed, the strain energy does 
not vary appreciably from that computed by direct averaging; 
therefore, even though the wing panel does exhibit the 
effects of bending and twisting, the internal strain energy 
computations do not seem to reflect this effect. The 
experimental analysis averages the inner and outer values 
of stresses before determining the panel strain energy; 
therefore, this method utilizes the midplane stresses that 
are the mean of the two surface stresses. It 1s significant 
to note that the stresses along the rosette legs must be 
tel lized £Or thie Computation. “The principal stresses 
cannot be used because the orientation of the inner and, 
outer stresses will not be parallel. 

The experimental result obtained utilizes different 
degrees of an averaging process. Obviously, this will not 
yield precisely the strain energy that actually exists in 


the panel; however, the partitioned panel strain energy 
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that was computed by summing the sub-panel energies must be 
as close as possible to the actual internal strain energy 
that exists. The further partitioning of the panel into 
smaller elements would require the installation of more 
strain gage rosettes. This would possibly further improve 
upon the accuracy of the experimental strain energy results; 
however, this approach was precluded by both time and 
monetary considerations. 

The optimum installation of the strain gage rosettes 
would have been to place them in the center of each of the 
fifteen sub-panel elements; however, the internal gages 
having been previously mounted and the desirability of 
mounting the inner and outer gages back to back precluded 
the optimum installation. Additionally, because of the 
lack of any stress concentrations, the stresses did not 
vary radically over the panel. For these reasons the 
rosette installations were considered adequate for this 


analysis. 
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CHAPTER VI 


COMPARISON OF RESULTS 


The comparison of the theoretical and experimental 
results, obviously, in the final analysis depends on how 
well the energy of the mathematical model approximates the 
true energy of the system. As stated in Reference 8, one 
Should not expect a correct result from either a mathema- 
tical model or a lumped averaging process similar to the 
empirical approach taken in this report unless the load 
distribution is uniform over the panel. The strain energy 
error becomes smaller if a more finely partitioned grid is 
used. Thus, when the panel is subdivided into smaller 
elements, the lLumped integral form of computing strain 
energy will increase in accuracy. 

The error occurring in the determination of the 


theoretical strain energy is obtained from the following: 


= percentage error (12) 


Utilizing the results obtained in these analyses, the 
percentage error in the theoretical strain energy computa- 
tion was found to be 2.94 per cent. 

It appears that there have been no gross errors made 
in the energy accounting procedures. The comparison made 
indicates a remarkably close agreement between experimental 


and theoretical results. 
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Superficially, one might expect the experimental strain 
energy to be greater than that computed from the idealized 
model. This would seem reasonable when one considers that 
the theoretical solution for a shear panel cannot possibly 
include the strain energies induced from bending, torsion, 
and any cross coupling effects that may exist in the actual 
thick milled skin panel; however, careful scrutiny of the 
idealized model reveals that the Poisson's effect is not 
included; and therefore, there is no mechanism to include 
the chordwise stresses. In fact, they are neglected. On 
the other hand, in the empirical analysis it is shown that 
appreciable chordwise stresses do exist. This is a contra- 
diction to the close agreement (2.94 per cent) in the 
comparison of the strain energies. 

From this discussion one must conclude that the 
stresses determined theoretically in Reference 3 must have 
been in excess of their true values. This excess of span- 
wise and shear stresses would compensate for the lack of 
chordwise stresses. This is the only logical explanation 
for the close comparison of the two analyses. 

In order to determine the validity of this reasoning, 
the principal stresses at station ¥,, = 34.50 inches were 
compared. The stresses in the model were: Ts =" Mocs: 
and 4. — S 596 psi. Similarly, the experimental stresses 


obtained from Figure 10 were: qd. = 525 psi and 


max 
~_ = Slo DSi. 
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As was predicted, the idealized normal stresses were 
considerably higher (26.30 per cent) than those obtained 
empirically. Similarly, the theoretical shear stresses were 
higher (13.25 per cent). In general, these overestimates of 
stress values compensated for the lack of chordwise stresses 
when the theoretical strain energy was computed. 

The location of the test panel on the wing did not lend 
itself to the direct application of severe bending and/or 
torsional loads. One can conclude, therefore, from the 
results of this report that when even thick skinned aircraft 
structures are analyzed using the conventional lLumped 
stringer and shear panel method, the analysis will be valid 
DEOVIded sehete areulossevere Dending Of torsion seoaeae 
applied directly to the panel. In general, the conventional 
method is valid, provided that the structure is sufficiently 
partitioned into a fine enough grid structure and the load 
distribution is not too irregular. 

In conclusion, one should not expect from an analysis 
based on the Lumped stringer and shear panel concept a 
satisfactory agreement with experiment for isolated load- 
ings, pabticularly when they are represented as point toreee 
in a theoretical analysis. In the panel studied and in" 
other similar cases there will always be ample structural 
material to further partition the model. This partitioning 
of a structure permits the desirable attenuation of the 


effect of load concentrations. These load concentrations 


can create discrepancies between computational and experi- 
mental results. 
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In an associated study, the wing center section is 
being included in the over-all wing analysis. The inclu- 
Sion of the energy contribution from the center section 
would not be expected to change the theoretical analysis in 
the vicinity of the particular panel considered in this 
report. This is because of the massive construction of the 
wing”in this specific” area. Even though the pivot rib ™is 
not infinitely rigid in the chordwise direction, it is 
nearly so and may be assumed to have zero flexibility. 
Reference 3 verifies this assumption by the close agreement 
between the theoretical and experimental results in the 


vicinity of the particular panel considered. 
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CHAPTER VII 


CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 


It can therefore be concluded that a valid comparison 
was achieved between theoretical and experimental strain 
energies of a thick milled skin aircraft structure. The 
accuracy of the analytical method is subject to the absence 
of concentrated loadings. Additionally, for a theoretical 
structure, a, fame clement ypart#tiening procedure is 
1 ecaimsede 

The analytical method provided remarkably accurate 
results (2.94 per cent error) and is a valid approach to 


aircraft structures of this specific type. 


Recommendations 


Le is recommended that G€ven though the results Of wri 
study are favorable, another analysis should be made on a 
panel that is subjected to more severe non=-planar loadings. 
Possibly a panel outboard of the intermediate rib could be 
used on the same wing. ‘The load would then be applied to a 
thinner section of the wing where there is less torsional 
rigidity, and the point of load application would be closer 
to the panel to be studied. Another possibility would be 
to obtain some means of applying single point loads to the 
same panel at the intersection of the intermediate rib and 


the spanwise beams. Still another method would be to apply 


52 


the same torque used in this analysis at the intermediate 
Yr ids’. 

An additional recommendation, that was previously dis- 
cussed in Chapter III, is that a study be conducted to 
determine the validity of a wing analysis that neglects the 


presence of any rib flanges. 


aS 
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APPENDIX A 
THEORETICAL DETERMINATION OF SHEAR FLOWS 
AT SECTIONS PERPENDICULAR TO CIB AT 


Vee — 98.7 and y= 74.3 (Reference 3) 


The theoretical. shear flows at sections perpendicular 
tet RewClLBaat Va 22. Ovondey = / 4-30 were calculated. in 
Order to determine the shear flows at both sections which 
would be used to produce the equivalent applied forces at 
the streamwise intermediate rib, Nhe = 81.98. The cross- 
sections are shown schematically in Figure Al. The thick- 
nesses were averaged desea values, and the web heights 
were taken as the mid-panel height between the upper and 
lower skins. 

A box structure with several cells will have one less 
redundant than the number of cells. In this case, there are 
five redundant webs. It is desired to write six equations 
in six unknown shear flows. This was done by equating the 
angle of twist of one cell with the remaining five, which 
gives five equations, and then writing an equilibrium of 
torsional moments equation. 

The first five equations were obtained by equating the 
angle of twist per unit length, oO , Of one cell with the 


other five using the well-kKnown expression for a box beam, 


9-5 ike . 
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FIG. Al 
SCHEMATIC DIAGRAM OF SECTIONS 
PERPENDICULAR TO THE C.LB. 
AT y, = 98.7 AND 743 


Upper surface 





Upper surface 





C13 


Note: dimensions in 
inches 
(ref. |O) 
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The equilibrium of torsional moments may be written as 


T= 2 2 Ang (A2) 


In using Equation A2 the summation is carried out 
SBOund scene .enctibe | perifiewm= Ot cach cell. The unit,.length, 
L, and the constant factor 2G drop out, leaving for cells a 


ang b;, 


The value of q for any exterior web of cell a is 
and for the interior ver is cee - Apt): Using the 


following abbreviations, 


Sen = 
Dbby camry te 
Sab L 


1} 


Equation A3 is rewritten, 
at 
te Sa ie Deb = — Cu Sub a 


tie terms Sao and Sop represent summations around the 


entire perimeter of their respective cells and Seb the 
value for the interior web. The enclosed areas were taken 


as the average web height times the distance between webs. 


Se 


Equating the angle of twist)per unit tengenot com 


to the remainirg cells yields the following equations: 


“e (Got bee = pot Sab) = 


Te Gue Suv - Got Sa Ger Sve) 
= 14 Gee ~ Got Ste - Gat Bcd) 
wale Bad ~ Get Bed - Get Sde) 
= Fe (Get Bec - git Bde -— 9% Bet) 
Ge bu - get Bef) 


Making appropriate substitutions, performing the indicated 
algebraic operations and rearranging, one obtains five 
equations in terms of the six shear flows with constant 
Geetficients. The required sixth equation was derived as 
previously described. 

The solution to these six Simultaneous equations was 
obtained from a digital computer using a FORTRAN program 
which is included in Reference 3. It utilizes Gauss's 
method of elimination with row pivoting and back substitu- 
tion and is designed to yield solutions for one or more 
column vectors forming the right hand side of the set of 


equations. This means that it would produce solutions for 


a 


one or more values of applied torque. The check values of 


q were hand calculated using an iterative procedure 
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suggested by Bruhn in Reference 9 and within slide rule 


accuracy agree well with the computer solutions. 
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Error- % of load 


Dynomometer reading Ib. x 1O"° 
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APPENDIX B 
CALIBRATION CURVES FOR DILLON DYNAMOMETERS 
AND HYDRAULIC PRESSURE SYSTEM 
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APPENDIX C 
CO-ORDINATE LOCATION OF STRAIN GAGES 


FS8U-3 WING PANEL 


Note: Co-ordinate origins are as follows: 


a O = Intersection of center section droop 
leading edge and the center line of 
the aircraft. 


a O = The aircraft center line. 


UPPER SKIN INSIDE GAGES 


Gage Number Gage Type 


83-84-85 
86-87-88 
89-90-91 
92-93-94 
95-96-97 
98-99-100 
Werk —TO 2-103 
10441 05—206 
107 -1L@8-109 
110-111-112 
113-114-115 
Jah6 —TAle/ =a ie 
ES) I 20a ak 
122-123-124 
25 =126= 127 
128-129-130 
40-41-42 





NOTES 


1. All strain gage rosettes werei mounted at the Naval 
Postgraduate School and oriented with one leg perpendicular 
to the CIB, the 45 degree gage pointing outboard and aft, 
and the third leg parallel to the CIB. 


2. All gages were manufactured by the Baldwin-Lima- 
Hamilton Corporation, taken from lot B-31 and rated at 
120.5 £ 0.5 ohms with a gage factor of 1.97 = 2 per cent. 


3. The low numbered rosette leg identifies the gage 
and is the leg parallel to the Load Reference Axis. 
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4. Gage No. 98 yields questionable results. 


UPPER SKIN OUTSIDE GAGES 


580-581-582 
583-584-585 
475-476-477 
586-587-588 
589-590-591 
992-593-594 
995-596-597 


I 96-399 —=600 
60L-—602 —cea 
604-605-606 
607-608-609 
610-611-612 
613-614-615 
BG — 6 leo 
619-620-621 
622-623-624 
625-626-627 





NOTES 


1. All strain gage rosettes were mounted at the Naval 
Postgraduate School and oriented with one leg perpendicular 
to the CIB, the 45 degree gage pointing outboard and aft, 
and the third leg parallel to the CIB. 


2. All gages were manufactured by the Baldwin-Lima- 
Ham lton Corporation, taken from lot B=31 and rated at 
120.5 + 0.5 ohms with a gage factor of 1.97 =© 2 per cent. 


3. The high numbered rosette leg identifies the gage 
and is the leg parallel to the Load Reference Axis. 
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APPENDIX D 


STRAIN GAGE INSTRUMENTATION 


The method of taking strain readings employed common 
Wheatstone Bridge circuitry as is shown in Figure Dl. The 
experimental test setup used the Budd/Datran Digital Strain 
Indicator (Model TC22). The Wheatstone Bridge and all 
associated circuits were incorporated within the indicator. 
The indicator interpreted the bridge unbalance and gave a 
visual display of the strain directly in units of micro- 
inches per inch. 

When the test runs were made, each of the twenty gages 
connected through the switching and balancing unit was 
zeroed on the digital strain indicator in the conventional 
manner prior to load application. At any time during the 
test run the calibration could be checked by merely pushing 
a butt@n on the counter. This gave a continual check and 
ensured the calibration of the counter at all times during 


the running of the test. 
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FIG. DI 
SCHEMATIC DIAGRAM 


STRAIN GAGE INSTRUMENTATION 






a 






V input 


PO 





D 
20 Channel 
Strain Gage 


Switching and 





Balance Unit 
Budd Co. 
Model C-IOT 


V input 


Digital Strain Indicator 
Budd Co. Model A-II0 


R, =R, = 120.5 ohms 
Re =Rz = 120.5 ohms 
Resistance tolerance: 5% 


Linearity tolerance: O5% 
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APPENDIX E 
EXPERIMENTAL PRINCIPAL STRESSES AND 


AXIS ORIENTATION 


The magnitudes of the principal stresses, maximum shear 
stresses and orientations were calculated for each rosette 
on both sides of the wing panel. The experimental strain 
gage readings are tabulated in Table EI. Calculations were 
Simplified through the use of the IBM 360 Digital Computer 
utilizing a FORTRAN IV program named "RECROSE" shown in 
Tables E2 and E3. The RECROSE results are listed in 
Table E4. 

The program was designed around the standard relation- 
ship between rectangular rosette readings and principal 
stresses found in any standard text (e.g., Wang in Refer- 


ence 7) on the subject. 





FiguserE. 


Finn E [698 4 oe Vee FR CHT 


os 
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rr... - Ei - &-e6F +e +h 
@ (Q--) se 


Tmax = 


ook (6-6) +[2é,-(E +6] 


Rel A 26.- (E,4 E.) 
Dp == avcta cae 


D. is the angle between leg number one and the maximum 
normal stress axis. A positive value indicates an angle in 
the direction of € 2: <a e and é. are the strain 
gage readings on their respective legs shown in Figure El. 
The number of the rosette leg parallel to the Load Reference 
Axis (CIB) is the identifying number for the rosette. 

The input to program RECROSE (Table EI) was designed to 
accommodate the above notation. Only the rosette identify- 
ing number is listed. This is also the gage number for 
which the first column of strains is listed. The second 
and third columns are then the values of the strain on the 
next two consecutive numbered gages comprising the rosette. 

This program was designed for 7079-T6 Aluminum with a 
Medulus of ElastYeity ‘of LOWS x 10° pounds per square inch 
ema a Poisson's ratio of 0.82. 

The results of RECROSE are listed in Table E4 and are 


identified by the rosette numbers which orient the™reference 


eels for Dz : 
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TABLE EI 
EXPERIMENTAL STRAIN READINGS 
(Input to RECROSE Program) 


336,000 sane lb. load 


INSIDE SKIN 


Rosette | | 4 i md 6 a. yh 
Number In, [C2 In |g mM. 
83 9 51 
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Computer 


Coded Name 


El 
E2 
E3 
SIGMAX 
SIGMIN 
SIGMA l 
SIGMA 3 
TAUMAX 
PHIPRR 
N 


\ 


E 


TABLE E2 


LIST OF SYMBOLS USED IN 


PROGRAM RECROSE 


Definition 
strain in identifying leg of rosette. 
strain in diagonal leg of rosette. 
strain in perpendicular leg of rosette. 
maximum principal stress. 
minimum principal stress. 
stress in identifying leg of rosette. 
stress in perpendicular leqg@ererosettre. 
maximum shearing stress. 


angle from Ey to Qr ae axis. 


The number of input rosettes. 


Modulus of Elasticity. 


KB , Poisson's ratio. 
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TABLE E3 


FORTRAN IV FROGRAM 


"RECROSE" 
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TABLE E4 
EXPERIMENTAL RESULTS 


(Output of Program RECROSE) 


S10) 6 eel) 
O47. O09 
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